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Abstract 

By considcring the decomposition of a genèric two qubit density matrix presented by Woot- 
ters [W. K. Wootters, Phys. Rev. Lett. 80 2245 (1998)], the robustness of entanglement for 
any mixed state of two qubit systems is obtained algebraically. It is shown that the robustness 
of entanglement is proportional to concurrence and in Bell decomposablc density matriccs it is 
equal to the concurrence. We also give an analytic expression for two separable states which 
wipe out all entanglement of these states. Since thus obtained robustness is function of the norm 
of the vectors in the decomposition we give an explicit parameterization for the decomposition. 

Keywords: Quantum entanglement, Robustness of entanglement, Concurrence 
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1 Introduction 

Quantum entanglement has recently been attracted much attention as a potential resource for 
communication and information processing [1, 2]. Entanglement is usually arise from quantum 
correlations between separated subsystems which can not be created by local actions on each 
subsystem. By definition, a mixed state p of a bipartite system is said to be separable (non 
entangled) if it can be written as a convex combination of product states 

i i 

where p^ and pf^ are states of subsystems 1 and 2, respectively. Although, in the case of pure 
states of bipartite systems it is easy to check whether a given state is, or is not entangled, the 
question is yet an open problem in the case of mixed states. Therefore having a measure to quantify 
entanglement of mixed states is likely to be valuable and several measures of entanglement have 
been proposed [3, 4, 5, 6]. 

One useful quantity introduced in [7] as a measure of entanglement is robustness of entangle- 
ment. It corresponds to the minimal amount of mixing with locally prepared states which washes 
out all entanglement. Authors in [7] have given analytical expression for pure states of binary 
systems. A geometrical interpretation of robustness is given in [8] and have pointed that two cor- 
responding separable states needed to wipe out all entanglement are necessarily on the boundary 
of separable set. Unfortunately, the above mentioned quantity as the most proposed measures of 
entanglement involves exteremization which are difncult to handel analytically. 

In this paper we evaluate robustness for a genèric two qubit density matrix. Our approach to 
finding robustness is based on the decomposition given by Wootters in [6]. Wootters in [6] has 
shown that for any two qubit density matrix there always exist the decomposition p = J2i \ x i) ( x i\ 
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in such a way that (xi\Xj) = Xióij where Aj are square roots of eigenvalues, in decreasing order, of 
the non-Hermitian matrix pp. We consider Wootters decomposition and show that for a genèric 
two qubit density matrix we can associate a tetrahedral with coordinates Pi = (xí\xí). It is shown 
that in terms of these coordinates the corresponding separable states form an irregular octahedral. 
Based on this we evaluate algebraically robustness of entanglement of a genèric two qubit density 
matrix by using a new norm defined by Wootters's basis. We also give an analytic expression 
for separable states that wipe out all entanglement and show that they are on the boundary of 
separable states as pointed out in [8]. Since in our approach of evaluation of robustness, one needs 
to know the norm of the vectors in the decomposition, we give an explicit parameterization for the 
Wootters decomposition to evaluate the norm explicitly. 

The paper is organized as follows. In section 2 using the Wootters's basis [6] we define a new 
norm to evaluate robustness. Robustness of entanglement of a genèric two qubit density matrix is 
evaluated in section 3. By giving an explicit parameterization for density matrix p, the norm of 
the vectors in Wootters decomposition of p is evaluated in section 4. The paper is ended with a 
brief conclusion in section 5. 

2 Physical norm based on Wootters's basis 

Here in this section we define a new norm to evaluate robustness of a genèric two qubit density 
matrix. To this aim we first review concurrence and Wootters's basis as presented by Wootters 
in [6]. Wootters in [6] has shown that for any two qubit density matrix p there always exist a 
decomposition 




(2-1) 
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called Wootters's basis, such that 

(xi\xj) = XiSij, (2-2) 
where Ai are square roots of eigenvalues, in decreasing order, of the non-Hermitian matrix pp and 

P = {(Jy ® <Jy)p*{(Jy <g> <Jy), (2-3) 

where p* is the complex conjugate of p when it is expressed in a Standard basis such as 
{|TT) > ltl)}> ÍIIT) > lli)} an d <7j/ represent Pauli matrix in local basis {|T)>|I)} • Based on this, 
the concurrence of the mixed state p is defined by max(0, Ai — À2 — A3 — À4) [6]. 

Now let us define states \xj} as 

\x' i } = ^=, for i = 1,2, 3, 4. (2-4) 

Then p can be expanded as 

p = J2*iWi)Wi\, (2-5) 

i 

and Eq. (2-2) takes the following form 

x i \x'^ = 8 i j. (2-6) 

Using the independency of Wootters's basis, one can expand an arbitrary 4x4 Hermitian matrix 
M in terms of them, that is, we can write 

M = J2a ij \x' i )(x' j \, (2-7) 

ij 

where the Hermiticity of M implies that a^- = a^. We now define norm of the matrix M by 



M\\ = ^\Tr(MM)\= \J2a^\, (2-í 
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where M is defined as Eq. (2-3). Obviously, the absolute value appearing in Eq. (2-8) guarantees 
the positivity of the above defined norm. In fact in cases that M is nonnegative matrix then M 
and MM is also nonnegative and absolute value sign in Eq. (2-8) can be neglected. 
Any Physical norm in quantum information should be at least invariant under local unitary trans- 
formation and below we prové that the norm defined in Eq. (2-8) is actually invariant under local 
unitary transformation. Suppose the matrix M is subjected to local unitary transformation defined 



as 



M — ► M' = (Ui <8> U 2 )M{U l U 2 ) ] = Y j ai j \x'l){x" j \, (2-9) 

x'l) defined as 



where \x") = {U\ <8> U 2 ) \x'j) are the new Wootters's basis with the corresponding 

x») = (a y a y ){U{ ® U* 2 ){a y ® a y ) i'S (2-10) 



which satisfy 

(x'lW'^=5 ih (2-11) 

where in Eq. (2-11) we have used the fact that (U T a y U)ij = —iekiUkiUij = — idet(U)eij = {a y )ij. 
Therefore using the definition of norm given in Eq. (2-8) and using Eq. (2-11) we have 



M = J|£4I ( 2 - 12 ) 



V % 3 

Now, comparing Eq. (2-12) with (2-8) we see that norm of M is invariant under local unitary 
transformations. 

With these considerations it is natural to define distance between two density matrices p\ and p 2 
as 



IIP1-P2II = yJ\Tr(( Pl - P2 )(pi ~ (2-13) 
which is invariant under local unitary transformations . 
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3 Robustness for two qubit density matrices 

According to [7] for a given entangled state p and separable state p s , a new density matrix p(s) 
can be constructed as, 

p( 8 ) = ——(p + 8 p a ), s >0, (3-14) 

S ~~T~ 1 

where it can be either entangled or separable. It was pointed that there always èxits the minimal s 
corresponding to p s such that p(s) is separable. This minimal s is called the robustness of p relative 
to p s , denoted by R{p \\ p s ). The absolute robustness of p is defined as the quantity, 

R(p || S) = min R(p \\ p s ). (3-15) 

Ps&S 

Du et al. in [8] have given a geometrical interpretation of robustness and pointed that if s in 
Eq. (3-14) is minimal among all separable states p s , i.e. s is the absolute robustness of p, then p s 
and p(s) in Eq. (3-14) are necessarily on the boundary of the separable states. 

Here in this section we obtain robustness for a genèric two qubit density matrix. Our method 
of evaluation of robustness is based on the decomposition of density matrix given by Wootters in 
[6]. By defining Pi = \Ki, where K- L = (x'Jx'j), then normalization condition of p leads to 

4 

Tr(p)=Y / Pi = h Pi>0. (3-16) 
i=i 

This means that with respect to coordinates Pi, the space of density matrices forms a tetrahedral. 

With respect to this representation separability condition Ai — À2 — A3 — À4 < takes the following 
form 

^_^_^_^< . (3-17) 

States that saturate inequality (3-17) form a plane called S\ (see Fig. 1). All states violating 
inequality (3-17) are entangled states for which Ai > À2 + A3 + A4. These states form an entangled 
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region with S\ as its separable boundary. There exist, however, three other entangled regions cor- 
responding with the cases that Xj is dominated (j = 2,3,4). These regions also define separable 
planes Sj. Four planes Si together with four planes S-, corresponding to Àj = 0, form an irregular 
octahedral corresponding to the separable states. This geometry is similar to that of Bell decom- 
posable states but here we have an irregular octahedral associated to separable states [9]. Figure 1 
shows a perspective of this geometry, where two separable planes S\ and S[ are shown explicitly. 

Now in order to obtain robustness of p suppose that a ray from p is drawn such that intersects 
the boundary planes of separable region at points p' s and p" s . Although p' s is necessarily on the 
plane Si, but p" is allowed to lie on any plane S[, S 2l S 3 or £4, where we evaluate robustness for 
each case separately. First consider the case that p" is lying on the plane S[. In this case p" can 
be written as a convex sum of three vertices of the plane 

p's = ^X" K)( x íl = a 2 cr 2 +a 3 cr 3 + 04 04, a 2 +a 3 + a 4 = l, (3-18) 

i 

where cjj are separable states that can be written as convex sum of two corresponding vertices of 
tetrahedral as 

" 3 - K^K t W> «I + K^Kl W> «I ' (3 - 20) 



and \'l are 



X'l = 0, (3-22) 



^7^ + 1^ <3 - M » 
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K = K^Kl + K^Kl- (3 " 25) 



By expanding p' s as convex sum of p and p" s 



P's = Y^- s (p + s Ps), (3-26) 



and also using the fact that the coordinates of p' s satisfy the equation 

P{ P> 2 í% Pi 



= 0, (3-27) 
K x K 2 K 3 K 4 ' V > 



after some àlgebra coordinates A^ of p' s = J2i K \ x i) ( x i\ can be written as 



( °-2 i «3 l «4 A \ 
1 (11 0.Q , íï/1 , (7 ' v ' 



K3+K4 "•" K 2 +if 4 2 



J «2 L 



y _ \^3+^4 X 2 +-ftT4 ^2+^3/ " ' ^ V£V2i-£V4 "-2 i-"- 3/ ,-g 90") 



)A 2 + f ( 



Cl'2 1 Q3 1 Q4 



0.2 1 Q3 1 Q4 



r)* + ¥( 



y _ \K 3 +Ka K2+K4, K2+K3 ) " z \j\3-t1\4 j\ 2 tA3y ,„ „„n 

3 a-) , 03 , a 4 , C ' \ > 



K3+K4 + K2+ÍÍ4 + K2+K3 



0-2 1 Q3 1 a 4 



>+¥( 





+ 


K 2 +K 3 


:+f 






ÍÏ2 


+ 




^3+^4 


K 2 +K 3 


; ^ 2 






012 


+ 


0-3 



y _ \K 3 +K4 ' K2+K4 ' K2+K3 J - ■ Z \J\3-tl\4 J\2+"4 y /o o 

4 q 2 1 03 , 04 _j_ Ç ' ^ ' 

X3+ÍÍ4 X2+ÍÍ4 ^ K2+K3 "T" 2 

where C = Ai — A2 — A3 — A4 is the concurrence of p. By using the above result and the definition 
of distance given in Eq. (2-13) one can evaluate robustness of p relative to p" s as 



, _ \\p-p's\\ Eifr-K) 2 Ç 

1 II n ' - n"\\ "Vv <\l - \"V 2a 2 , 2a 3 , 2a 4 ' V° 

Wfs l·'sW V "i) K3+K4 K2+K4 ~*~ K2+K3 

Analogue to the above method one can evaluate robustness of p for the case that p" s lies on the 
plane S2. In this case p" s can be expanded as convex sum of three vertices of the plane 

pg = &i<ti + 630-3 + Ò40-4, bi + 63 + 64 = 1, (3-33) 

where 

" = k^tkí w> « 1 + ifrbï l4> ^ • (3 - 34) 
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and 03 and 04 are defined in Eqs. (3-20) and (3-21). Then after some àlgebra we obtain the 
corresponding robustness as 

c 

s 2 = 2b, | 2b 4 - ■ ( 3 " 35 ) 

K 2 +K 4 "T K 2 +K 3 

Similarly in cases that separable state p" are on the planes S3 and 1S4 we obtain relative robustness 
of p as 

5 3 = 2c 2 C 2c a ~ > ( 3 " 36 ) 

K3+K4 ^ K2+K4 

and 

54 = 2d 3 | 2d4~ ' ( 3 " 37 ) 

K 2 +K 4 K 2 +K 3 

respectively. Equations. (3-32), (3-35), (3-36) and (3-37) show that in order to achieve the minimum 
robustness it is enough to consider the case that separable state p" lies on the plane S[ . With this 
consideration we are now allowed to choose coefhcients in such a way that Eq. (3-32) becomes 
minimum. It is easy to see that this happens as long as the coefhcient au corresponding to the term 
min(iíj + Kj) becomes one. Therefore robustness of p relative to p" s is 

s = mm( ^ + ^ g (3-38) 

which is the main result this work. Here the minimum is taken over all combination of Ki + Kj 
for i,j = 2,3,4. Equation (3-38) implies that for two qubit systems robustness is proportional to 
the concurrence. We see that the minimum robustness given in Eq. (2-3) corresponds to aj = 5^-, 
therefore, by using Eq. (3-18) we get the following result for p" s 

Ps = °k- (3-39) 

Also by using = 5ik in Eqs. (3-28) to (3-31) one can easily obtain the coordinates of separable 
state p' s . 
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As we will show in the next section the Bell decomposable states correspond to the K i = 1 for 
i = 1,2,3,4, therefore in Bell decomposable states Eq. (3-38) implies that the robustness is equal 
to the concurrence. 

Now we have to show that thus obtained robustness is minimum over all separable states. Vidal 
et al. in [7] have shown that R(p \\ p s ) is a convex function of p s . This means that any local 
minimum is also the absolute one, thus in order to find the absolute minimum of R(p \\ p s ) as a 
function of p s it is enough to find local minimum [7]. As the above mentioned argument shows 
robustness given in Eq. (3-38) is minimal relative to all separable states of the tetrahedral. In the 
rest of this section we want to show that it is indeed minimum over all separable states. To this 
aim let us consider following pseudomixture for the density matrix p given in Eq. (2-5) 

p=(l + s)p' s -sp" s , (3-40) 

where p' s and p" s are two separable states with following decomposition 

p's = E A ÍK><^I + E a iM > < 4 1 > (3- 41 ) 
p'; = E A "K)(^l + E^l^)(4l' ( 3 - 42 ) 

i i,j 

where A and B are, respectively, off-diagonal Hermitian matrices A = (a^) and B = (bij). It 
follows from Eqs. (3-40) to (3-42) that the following equations should hold 

A i = (l + s)AÍ-sAÍ', (3-43) 

(1 + s) aij - s bij = 0. (3-44) 
Now using Eqs. (2-13), (3-41) and (3-42) robustness of p relative to p" s is defined by 



II P ~ Ps 

s = 



p' s -pU\ \ 



Tr(p - p' s ){p - p' s ) 



Tr(p< s -p<l)(p< s -pi>) 
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s = 



Zi(\i-\d 2 + Tr(AA*) 



\| E t (K-^) 2 + Tr(A~B)(A-By 



Zi(Xi-W + Tr(AA*) 



(3-45) 



where in the last line Eq. (3-44) have been used. in order to show that thus obtained robustness 
is minimum over all separable states we need to show that it is local minimum [7]. To do so, it is 
straightforward to see that the robustness s given in Eq. (3-45) reduces to 



which takes minimum vàlues given in Eq. (3-38) as the earlier argument of this section indicates. 
Since as long as both the numerator and denominator under the radicals have the same sign, s 
would be independent of the perturbation of corresponding separable states. While in cases that 
numerator and denominator posses different sign then s would be perturbation dependent. It is 
trivial to see that the numerator and denominator should change sign simultaneously, hence s 
would be independent of perturbation. Since both the numerator and denominator of robustness 
s given in Eq. (3-45) are continuous function of the matrix elements of matrices A and B, so we 
expect that s is also continuous function of them. Now suppose that separable states p' s and p" 
given in Eqs. (3-41) and (3-42) are obtained by small perturbation of corresponding separable 
states on planes S\ and S[. This means that matrix elements of A and B are infinitesimal , i.e., 
Tr(AA*) is infinitesimal. Though Tr(AA*) can be negative but both numerator and denominator 
of Eq. (3-45) remain positive. Therefore for small enough perturbations the robustness s is equal 
to (3-46) and it is independent of perturbation. Hence if the numerator changes the sign under the 
appropriate continuous perturbation while the denominator remains positive, obviously robustness 
s would vanishes at perturbation corresponding to the change of sign of the numerator. Similarly if 




(3-46) 
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the denominator changes sign while the numerator remain positive, the robustness becomes infinite 
at the corresponding perturbation. In either case the robustness s as a continuous function of 
parameters of perturbation will jump discontinuously to zero or infinite which is not possible. 

In summary as we see the continuous perturbation of separable states which minimum robust- 
ness can not affect the robustness. This means that the obtained robustness given in Eq. (3-38) is 
local minimum thus according to Ref. [7] it is global minimum. 



4 Evaluation of Ki via the explícit parameterization of the density 
mat r ix 



In the previous section we evaluated robustness for a genèric two qubit density matrix. As we see our 
approach of evaluation of robustness is based on the decomposition (2-1). As we see thus evaluated 
robustness is function of concurrence and also JQ, i.e., norm of vectors of the decomposition. Here 
in this section we introduce an explicit parameterization for \xí), in order to evaluate K{. To this 
aim we define matrix X and X as 



x = (|xi) , |4) , 14) > \ x 'a)) > 



and 



X 











4), 


4) ) 


4) ' 





(4-47) 



(4-48) 



respectively. Therefore Eq. (2-2) takes the following form 



X ] X = X T a y ® G y X = I. 



Symmetric matrix a y ® o y can be diagonalized as 



(4-49) 



(4-50) 
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where O is an orthogonal matrix defined by 



and i] is the diagonal matrix 



V 



( \ 
10 1 

11 

1-10 



\ 1 ~ l J 



( \ 

% 

10 

i 



1 

Using Eq. (4-50) we can rewrite Eq. (4-49) as 



V 



(4-51) 



(4-52) 



Y T Y = L 



(4-53) 



where 



Y = r] OX. (4-54) 

Equation (4-53) shows that Y is a complex 4-dimensional orthogonal matrix. This means that a 
given density matrix p with corresponding set of positive numbers Àj and Wootters's basis can 
transforms under SO (4, c) into a genèric 2x2 density matrix with the same set of positive numbers 
but with new Wootters's basis. This implies that the space of two qubit density matrices can be 
characterize with 12-dimensional (as real manifold) space of complex orthogonal group SO(4,c) 
together with four positive numbers \. Of course the normalization condition reduces number of 
parameters to 15. 
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As far as entanglement is concerned the states p and p' are equivalent if they are on the same 
orbit of the group of local transformation, that is, if there exist local unitary transformation U\ <8> U2 
such that p' = {U\®U2)p(Ui®U2)\ where U\ and U2 are unitary transformations acting on Hilbert 
spaces of partides A and B, respectively. 

It can be easily seen that under the above mentioned local unitary transformations of density 
matrix p, the matrix X transforms as 

X -> X' = (Ui ® U 2 )X. (4-55) 

It is worth to mention that X' also satisfy Eq. (4-49). To show that this is indeed the case, we 
need to note that X' 7 u y <8> o y X' = X T {U^ a y U\) <8> (C/J \r y U2)X . By using (a y )ij = —itij we get 
(U T a y U)ij = —ieuUkiUij = —idet(U)eij = {cr y )ij, where last equality comes from the fact that the 
unitary matrix belongs to SU(2). This implies that 

X' T Gy <g> GyX' = I. (4-56) 

By defining Y' as 

Y' = r]OX', (4-57) 
one can easily show that Y' is also satisfies orthogonality condition 

Y' T Y' = I. (4-58) 

Now by using Eq. (4-57) and inverting Eq. (4-55), we can express Y' in terms of Y 

Y' = (ï?0)(C/i ® U 2 ){r]0)- l Y. (4-59) 

Now by using the fact that (r)0)exp(Ui ®U.2){r}0)~ 1 = exp((ry 0)(U\ ® Z^Xrç O) -1 ) and using 
the explicit form for generators (U\ <8> W2) of local group, one can after some algebraic calculations 
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see that (rjO)(Ui ®U2){j}0)~ 1 is real antisymmetric matrix. This means that under local unitary 
transformations matrix Y transforms with 50(4, r) group. So we can parameterize the space of 
two qubit density matrices as 6-dimensional coset space 50(4, c)/50(4, r) together with 4 positive 
numbers Àj, which again normailzation condition reduces the number of parameters to 9. 

Below in the rest of this section we will obtain an explicit parameterization for a genèric two 
qubit density matrix. First note that we can decompose coset 50(4, c)/50(4, r) as [10] 



50(4,c)/50(4,r) 



fSO(2,c) SO(2,c) 



SO{A,c) 

50(4,r) ~ 50(2, c)/SO{2, r) <g> 50(2, c)/SO(2, r) \SO{2,r) " SO(2,r) J ' 



(4-60) 



that is, coset representative Y can be decomposed as Y = Y\Í2- One can easily show that coset 
representative of SO(2,c)/SO(2,r) has the following form 



exp 



o i4 

-i<f> 



cosh (f> i sinh <p 
— isinh0 cosh</> 



(4-61) 



Thus Y2 can be written as 



Yo 



On the other hand Y\ can be evaluated as 



cosh (f>i i sinh <j)\ 









—i sinh 0i cosh 4>i 






cosh 4>2 i sinh (f>2 









—i sinh 4> 2 cosh 2 



(4-62) 



Y\ = exp 






iB 


-ÍB T 






cosh VBB T 



• sinh VB T B _gT 



iB 



sinhVB T B 



cosh \/B T B 





\ 

iC 


Vi + cc T 


-iC T 




Vi + c T c J 



(4-63) 



where B is a 2 x 2 matrix and in the last step we used C = B smh , == ^ . Now using the singular 

v5 B 
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value decomposition C = O1DO2 , Eq. (4-63) becomes 



Y-i = 



o 1 VTTd^oJ 



V 



iOxDOl 



\ 



2 Vl + D 2 2 T j 



(4-64) 



-i0 2 DOf 

where D is a non-negative diagonal matrix. It can be easily seen that Eq. (4-64) can be decomposed 
as 



Yi = 



( 


°1 


( l 




Oi 


VI + D 2 


iD 


{ 


o 2 J 




Vl + D 2 



By combining Eqs. (4-62) and (4-65) we get 
Y = 



\ í 



1 


] 


( , 




Oi 




iD 


{ 


o 2 ) 




VI + D 2 



Finally using parameterization given in Eq. (4-62) we get 



V 



oi 








ol 


Oi 








o 2 



(4-65) 



(4-66) 



Y = 



( 


cosh 61 i sinh 0\ 





\ 


/ 


cosh Çi 





i sinh Çi 







—i sinh 9\ cosh 0\ 









cosh Ç2 





i sinh Ç2 







cosh 62 i sinh 62 






— isinh £1 





cosh£i 





V 


—i sinh 62 cosh (9 2 


/ 


V 





-i sin Ç 2 





cosh Ç2 



V 



cosh (f>\ i sinh 0i 









— ísinh0i cosh 0i 






cosh (p2 i sinh 02 









—i sinh 02 cosh 02 



(4-67) 



Using above results and Eq. (4-47) and (4-54) we can evaluate the states \xí) as 
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'A_2 



'A_3 



^ — (sinh £1 sinh 9 2 cosh 0i + sinh £ 2 cosh 9 2 sinh <j>i) — i (cosh £1 cosh 9\ cosh 0i + cosh £ 2 sinh 0i sinh 0i ) 

— (sinh £1 cosh 9 2 cosh 0i + sinh £2 sinh 9 2 sinh 0i ) — i (cosh £1 sinh 9\ cosh 0i + cosh Ç2 cosh #1 sinh 0i ) 

(sinh £1 cosh 9 2 cosh 0i + sinh £2 sinh 9 2 sinh 0i ) — i (cosh £1 sinh 9\ cosh 0i + cosh £2 cosh d\ sinh 0i ) 

(sinh £1 sinh 9 2 cosh <p\ + sinh £ 2 cosh 9 2 sinh 0i ) — i (cosh £1 cosh 9\ cosh 0i + cosh £ 2 sinh 9\ sinh 0i ) 

(4-68) 

'cosh ^ cosh 6*i sinh </>! + cosh £ 2 sinh 6*i cosh — z (sinh £1 sinh 9 2 sinh 0i + sinh £ 2 cosh 2 cosh X ) 

cosh^i sinh í^i sinh0i + cosh £2 cosh#i cosh0i) — i (sinh £1 cosh #2 sinh0i + sinh £2 sinh #2 cosh0i) 

cosh £1 sinh 9\ sinh 0i + cosh £2 cosh #1 cosh 0i ) + i (sinh £ 1 cosh 9 2 sinh 0i + sinh £ 2 sinh #2 cosh 0i ) 

cosh^i cosh í?i sinh 0i + cosh £2 sinh#i cosh0i) + i(sinh£i sinh $2 sinh0i + sinh £2 cosh #2 cosh0i) j 

(4-69) 

sinh Çi cosh 9\ cosh (f> 2 + sinh £2 sinh 9\ sinh 02 ) — i (cosh £1 sinh 9 2 cosh 02 + cosh £ 2 cosh #2 sinh 2 ) 

sinh Çi sinh #1 cosh <p 2 + sinh £ 2 cosh 9\ sinh 02 ) — i (cosh £1 cosh # 2 cosh 02 + cosh £ 2 sinh 9 2 sinh 2 ) 

sinh Çi sinh #1 cosh 2 + sinh £ 2 cosh #1 sinh 2 ) + i (cosh Çi cosh 9 2 cosh 2 + cosh £ 2 sinh # 2 sinh 2 ) 

y (sinh £1 cosh #1 cosh 2 + sinh £ 2 sinh #1 sinh 2 ) + i (cosh £1 sinh 9 2 cosh 2 + cosh £ 2 cosh # 2 sinh 2 ) 

(4-70) 

(cosh £1 sinh # 2 sinh 2 + cosh £ 2 cosh # 2 cosh 2 ) + i (sinh £1 cosh #1 sinh 2 + sinh £ 2 sinh 9\ cosh 2 ) 
(cosh^i cosh # 2 sinh 2 + cosh £2 sinh #2 cosh0 2 ) + i (sinh £1 sinh#i sinh0 2 + sinh £2 cosh#i cosh 02 ) 

— (cosh £1 cosh 9 2 sinh 2 + cosh £ 2 sinh 9 2 cosh 2 ) + i (sinh £ 1 sinh 9\ sinh 2 + sinh £ 2 cosh 9\ cosh 2 ) 



— (cosh £1 sinh 9 2 sinh 2 + cosh £ 2 cosh 9 2 cosh 2 ) + i (sinh £1 cosh 9\ sinh 2 + sinh £ 2 sinh #1 cosh 2 ) 

(4-71) 



Equations (4-68) to (4-71) together with normalization condition J2Í=i ( x i\ x i) = 1 gi ye a parame- 
terization for a genèric orbit of two qubit density matrix up to local unitary group. We are now 
in position to evaluate K{. By using Eqs. (4-68) to (4-71) after some àlgebra we arrive at the 
following results 

K\ = cosh 262 (sinh 2 £1 cosh 2 0i + sinh 2 Ç2 sinh 2 0i) 
+ cosh 26*i (cosh 2 Çi cosh 2 <f>\ + cosh 2 £2 sinh 2 0i) 
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+ sinh 20i (sinh Çi sinh £2 sinh 202 + cosh £1 cosh £2 sinh2#i), (4-72) 

K 2 = cosh 26> 2 (sinh 2 £1 sinh 2 fa + sinh 2 £ 2 cosh 2 fa) 
+ cosh 26*i (cosh 2 £1 sinh 2 fa + cosh 2 £2 cosh 2 0i) 

+ sinh 2 fa (sinh £1 sinh £2 sinh 29 2 + cosh £1 cosh £2 sinh 29 1), (4-73) 

K3 = cosh 2^i (sinh 2 £1 cosh 2 02 + sinh 2 £2 sinh 2 2 ) 
+ cosh 2^2 (cosh 2 Çi cosh 2 2 + cosh 2 £2 sinh 2 02 ) 

+ sinh 202 (sinh Çi sinh Ç2 sinh 2^i + cosh £1 cosh £2 sinh 26*2), (4-74) 

K 4 = cosh 29 1 (sinh 2 £1 sinh 2 2 + sinh 2 £ 2 cosh 2 2 ) 
+ cosh 2^2 (cosh 2 £1 sinh 2 02 + cosh 2 £2 cosh 2 02) 

+ sinh 202 (sinh Çi sinh Ç2 sinh 2^i + cosh £1 cosh £2 sinh 26*2)- (4-75) 

As an example let us consider Bell decomposable states p = 2~2í=iPí I Y^) (V'il) where > 
0) Si K = 1- For these states by choosing 9\ = #2 = £1 = £2 = 01 = 02 = we get Aj = and 
states \xí) are given by 





fPÏ\fa) , 


hM = 


^(ITT) + IU)), 


(4-76) 


M = y 


/P2\tp2) , 


hk> = 


^(ITD + IIT)), 


(4-77) 


X3) = ~i\ 


/P3 1^3) , 


l^ 3 ) = 


^(ITD-IIT)), 


(4-78) 


|x 4 ) = ■ 


s/PA |V>4> , 


hM = 


^(ITT)-IU)). 


(4-79) 



So it can be easily seen that for these states we get = 1 for i = 1, 2, 3,4. 

5 Conclusion 

In this work we have obtained robustness of entanglement of a genèric two qubit density matrix. 
Our approach to obtain robustness is based on the decomposition of a genèric two qubit density 
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matrix presented by Wootters. We have shown that the robustness of entanglement is proportional 
to concurrence and in Bell decomposable density matrices it is equal to the concurrence. We also 
present an analytical expression for two separable states that wipe out all entanglement of these 
states. We show that robustness is function of the norm of the vectors of the decomposition so we 
present an explicit parameterization for the decomposition. 
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Figure 1: The space of a genèric two qubit density matrix is represented by a tetrahedral. 
Vertices P{ for i = 1,2,3,4 are correspond to pure states defined by p = Aj|a^)(a^|. Irregular 
octahedral corresponds to separable states. Separable planes S\ and S[ are shown explicitly. 
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